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Abstract 

We have calculated the next-to-leading order perturbative QCD correc- 
tions to the photon energy spectrum in radiative T decays. The higher-order 
corrections significantly modify the shape of the spectrum, in particular at 
large photon energies, and thereby reduce the discrepancy between experi- 
mental data and previous leading-order predictions. The next-to-leading order 
calculation of the photon energy spectrum allows a more reliable determina- 
tion of the strong coupling constant from radiative T decays by restricting 
the analysis to the region of the energy spectrum that can be described by 
NLO perturbation theory. 
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1. Introduction. The calculations of heavy quarkonium decay rates are among 
the earliest applications of perturbative QCD they have been used to extract 
information on the QCD coupling at scales of the order of the heavy-quark mass. A 
consistent and rigorous framework for treating inclusive quarkonium decays has re- 
cently been developed, superseding the earlier non- relativist ic potential model. The 
factorization approach p[ is based on the use of non-relativistic QCD (NRQCD) 
to separate the short-distance physics of heavy-quark annihilation from the long- 
distance physics of bound-state dynamics. The quarkonium decay rate can be ex- 
pressed as a sum of terms, each of which factors into a short- distance coefficient 
and a long-distance matrix element. The short-distance coefficients are determined 
by the annihilation cross section for an on-shell QQ[n\ pair in a colour, spin and 
angular-momentum configuration n, and can be calculated perturbatively in the 
strong coupling as{mQ). The probability to find a QQ pair in the quarkonium at 
the same point and in the configuration n is parametrized by a non-perturbative 
matrix element (C^), which is subject to the power counting rules of NRQCD and 
which can be evaluated numerically using NRQCD lattice simulations Q]. It is thus, 
in principle, possible to calculate the annihilation decay rates of heavy quarkonium 
from first principles, the only inputs being the heavy-quark mass uiq and the QCD 
coupling asimo). Fixing the heavy-quark mass from, for example, sum rule cal- 
culations for quarkonia, or considering ratios of different decay channels where the 
quark-mass dependence and the long-distance factor cancel, the analysis of quarko- 
nium decay rates can provide information on the QCD coupling at the heavy-quark 
mass scale. 

In this letter, we focus on the reaction T — > 7 -|- X, which is used for as determi- 
nations from the quarkonium system and which provides detailed information on the 
decay dynamics through the analysis of the photon energy spectrum. We present 
the first calculation of the next-to-leading order perturbative QCD corrections to 
the photon spectrum. Besides being an interesting problem of perturbative QCD 
per se, the inclusion of higher-order corrections allows a NLO determination of the 
strong coupling constant, avoiding those regions of the energy spectrum that are 
affected by non-perturbative contributions.^] 

2. The leading-order calculation. At leading order in the non-relativistic 
velocity expansion in f ^, radiative T decays proceed through the annihilation of a 
colour-sing let n = ^S{^^ bb pair (in spectroscopic notation) .0 The photon energy 

-"^In the absence of a NLO prediction for the photon spectrum, model approaches have been 
used to extrapolate the experimental data into the region of low photon energies to obtain the 
total decay rate, introducing uncontrolled theoretical uncertainties. The model dependence of the 
corresponding NLO determination of as can be avoided by restricting the analysis to the region of 
the energy spectrum that is described by NLO perturbation theory. 

^At leading order in the velocity expansion, the NRQCD description of 5'-wave quarkonium 
decays is equivalent to the non-relativistic potential model, where the non-perturbative dynamics is 
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spectrum receives contributions from two different mechanisms. At leading oder, 
one has 



, (1) 



where Q < x, = / mi) < 1 . The first term on the right-hand side denotes the direct 
cross section, where the photon is radiated off a heavy quark, Fig. |l](a), while the 
second term represents the contribution from gluon fragmentation into a photon. 
Fig. 0(b). The fragmentation contribution, although of 0{aat) in the coupling 
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Figure 1: Generic leading-order Feynman diagrams for radiative T decay: (a) direct 
contribution, (b) fragmentation contribution. 

constant, is enhanced by a double logarithmic singularity ~ ln^(m^/A^) ~ l/Og 
arising from the phase-space region where both the light-quark-antiquark splitting 
as well as the photon emission become coUinear 0]. The fragmentation function 
is sensitive to non-perturbative effects and has to be extracted from experiment. A 
comparison between the direct component and the fragmentation contribution to 
the photon spectrum, performed at leading order, reveals that the fragmentation 
process is important in the low-x-y region, but suppressed with respect to the direct 
cross section for x, ^ 0.3 It is, however, important to point out that the total 

decay rate, including fragmentation processes, r(T 7X) = /(} dx^dF/dx^, is not 
an infrared-safe observable. For — 0, the fragmentation contribution rises like 
characteristic of the soft bremsstrahlung spectrum, and cannot be integrated 
down to x-y = 0. Therefore, in perturbation theory, only the photon energy spectrum 
can be calculated for x^ 7^ 0, after coUinear singularities have been factorized into 
the fragmentation functions. 

In the following, we will focus on the direct contribution of the energy spec- 
trum, which dominates in the region x, ^ 0.3 where accurate experimental data are 

described by a single long-distance factor related to the bound-state's wave function at the origin. 
In the case of P-wave quarkonia, the potential-model calculations at 0{a^) encounter infrared 
divergences, which cannot be factored into a single non-perturbative quantity. Within the NRQCD 
approach, this problem finds its natural solution since the infrared singularities are factored into 
a long-distance matrix element related to higher Fock-state components of the quarkonium wave 
function. 

■^At next-to-leading order, the fragmentation contribution involves also quark fragmentation, 
which is much harder than gluon fragmentation and may be important at larger values of x-y. 
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available. At leading order, one finds 

^ dr^g(T 7X) ^ 2 / _ . _ 2 (1 - ^7)^ 1,(1 _ , ) 

I ^ ~ -^7 _|_ 3^7(1 ~ 3^7) I 

To very good accuracy, Eq. (^) can be approximated by a linear spectrum: l/rjj° 
dFj^/dx^ ~ 2x^. The corresponding leading-order decay rate is given by 

r^°(T - 7X) = ^ a el {n^ - 9) ^^i^, (3) 

where e?, is the charge of the b quark. Up to corrections of C(f the non-perturbative 
NRQCD matrix element is related to the T wave function at the origin through 
{Oji^Si)) fti 2Nc I (^(0) p, with Nc = 3 the number of colours. 

The linear rise with predicted by LO perturbative QCD, Eq. (|^), is not sup- 
ported by the experimental data pf. In particular for ^ 0.7 the discrepancy 
is significant. The leading-order estimate of the photon energy spectrum demands 
several theoretical refinements: the inclusion of (i) higher-order perturbative QCD 
corrections; (ii) relativistic corrections due to the motion of the b quarks in the T 
bound state; and (iii) hadronization effects that may become important close to the 
phase-space boundary. The relativistic corrections to the photon energy spectrum, 
including colour-octet contributions, have been analysed in Although poten- 

tially sizeable, they are concentrated in the upper and lower end-point region and 
do not significantly modify the shape of the spectrum in the intermediate energy 
range 0.4 <^ x^ <^ 0.9. Hadronization effects have been studied in the context of 
different models and may become important when ^ 1. However, these 



non-perturbative contributions are not expected to modify the energy distribution 
significantly in the intermediate region x-y ~ 0.7. The calculation of the next-to- 
leading order perturbative QCD corrections to the photon spectrum is the subject 
of this letter. 

3. Next-to-leading order QCD corrections. Generic diagrams that build 
up the decay rate in 0{aa^) are drawn in Fig. ^ The evaluation of the next-to- 
leading order corrections involves two parts, the virtual corrections, generated by 
virtual particle exchange, and the real corrections, which originate from real gluon 



radiation as well as from gluon splitting into light-quark-antiquark pairs [jTT . 

The virtual amplitude consist of 105 diagrams, including self-energy and vertex 
corrections for photon and gluons (2a,b) as well as a large number of box diagrams 
(2c-e). The evaluation of these amplitudes has been performed in the Feynman 
gauge. The ultraviolet and infrared singularities have been extracted using dimen- 
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Figure 2: Generic next-to-leading order Feynman diagrams for direct radiative T decay. 

sional regularization, leading to poles in e = (n — 4)/2. The Feynman integrals con- 
taining loop momenta in the numerator have been reduced to a set of scalar integrals 
using an adapted version of the reduction program outlined in |T^. This program 
has been extended to treat n-dimensional tensor integrals with linear-dependent 
propagators to account for the special kinematical situation in the non-relativistic 
approximation to T decays. Using the Feynman parameter technique, we have ob- 
tained analytical results for the scalar integrals; these will be presented in a forth- 
coming publication. The masses of the riif = 4 light quarks have been neglected, 
while the mass parameter of the h quark has been defined on-shell. We have carried 
out the renormalization of the QCD coupling in the MS scheme, including nif = 4 
light flavours in the corresponding j3 function. The bottom quark has been decou- 
pled from the evolution of by subtracting its contribution at vanishing momentum 
transfer [jl3| . The exchange of Coulombic gluons in diagram (2e) leads to a singular- 
ity ~ tt"^ /2f3R, which can be isolated by introducing a small relative quark velocity 
I3r. The Coulomb-singular part can be associated with the interquark potential of 
the bound state and has to be factored into the non-perturbative NRQCD matrix 
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element. Only the exchange of transverse gluons contributes to the next-to- leading 
order expression for the short-distance annihilation rate. 

The real corrections are generated by gluon radiation off quarks and gluons 
(2f,g) and by gluon splitting into light- quark-ant iquark pairs (2h). The 54 diagrams 
have been calculated in the Feynman gauge. For the sum of the gluon polarization 
we have used Y,^^l^u = —Qnv and the unphysical longitudinal gluon polarizations 
have been removed by adding ghost contributions. The real emission cross section 
contains singularities in the limit where the gluon splitting into gluon or light-quark- 
antiquark pairs becomes soft or collinear. We have used the subtraction method 
|T^ as formulated in to extract the singular parts of the real cross section in 



n dimensions. The method of [|15| is based on the fact that the soft and collinear 
limits of a real emission matrix element can be expressed in terms of a convolution 
of process- independent splitting functions and the leading-order matrix element. 
Subtracting from the real matrix element counterterms constructed accordingly, 
the infrared and collinear singularities cancel and the phase-space integration can 
be performed numerically in four dimensions. The counterterms can be integrated 
analytically in n dimensions over the phase space of the extra emitted parton, leading 
to poles in e = (ra — 4)/2. When the integral of the subtraction counterterm and the 
virtual cross section are combined, the infrared and collinear singularities cancel. 
The analytical result for the matrix element squared has been implemented in a 
Monte Carlo integration program, so that not only the inclusive decay rate but any 
distribution can be generated at next-to-leading order. 

As mentioned above, owing to the presence of fragmentation processes at small 
photon energies, the total decay rate is not an infrared-safe observable. We will 
nevertheless present a next-to-leading order result for the perturbatively calculable 
direct contribution to the total decay rate, in order to compare our calculation with 
a previous NLO estimate of that quantity. We find 



rS^°(T 7^) = 1 - ^ (l.71 + /5o(n,) In 



(4) 



where /5o(?T.if) = (llA^^c — 2nif)/3. The theoretical uncertainty due to the numerical 
phase-space integration has been reduced to ^ 0.5%. The new, more accurate result 
for the total rate is consistent with the previous calculation |T^, which yielded 
Fji^^/rjj'^ = (1 — Q;s(1.67 ± 0.36)/7r) for /i = 2m5, with a considerable theoretical 
uncertainty of ~ ±20% coming from the numerical evaluation of the loop integrals.^ 
The central part of this work is the calculation of the photon energy spectrum 
at next-to-leading order. The result is presented in Fig. |^. The QCD corrections 
significantly flatten and deplete the spectrum for ^ 0.75. The shape of the 
NLO spectrum is in qualitative agreement with the experimental data and indicates 

prescription for fixing the renormalization scale /i has been given in [T^ . The analysis 
implies that it is advantageous to consider ratios like r(T — > ^gg)/T{T — > ggg) for as extractions. 
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Figure 3: Photon energy spectrum for direct radiative T decays in leading and next-to- 
leading order (og = 0.2). 



that the discrepancy between theory and experiment at large photon energies can 
be reduced by the inclusion of higher-order QCD corrections. A proper comparison 
with the recent measurements requires the inclusion of experimental aspects, such as 
energy resolution and efficiency, and will be presented in a forthcoming publication. 

Finally, let us briefly comment on the general structure of the higher-order cor- 
rections near the upper end-point of the photon energy spectrum. It has been argued 
that potentially large logarithms ln(l — Xy) associated with the imperfect cancel- 
lation between real and virtual emission of soft gluons as — > 1 may contribute 
to all orders in perturbation theory [|TB|. The resummation of these soft-gluon ef- 

x^)) near the 
a 



fects may then give rise to a Sudakov suppression ~ 1/ exp(as In (1 
end-point of the photon energy spectrum. In contrast to the previous claim 
recent analysis [|l9l flnds that the logarithms of (1 — x-y) cancel at each order in the 
perturbative expansion, as long as the T decay proceeds through the annihilation 
of a colour-singlet bb pair. While the numerical NLO calculation of the spectrum as 
presented in Fig. ^ is not suited to deflnitely resolving that issue, it may be possible 
to derive, from the full NLO result, a simple analytic expression of the NLO cross 
section near the end-point region. We leave such a calculation for a forthcoming 
publication. 
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In conclusion, the photon spectrum in radiative T decays is a very interesting 
laboratory to study the effects of higher-order perturbative QCD corrections. The 
next-to-leading order calculation presented here shows that the intermediate region 
of the energy spectrum can be adequately accounted for by perturbative QCD. 
Fragmentation processes, relativistic corrections and hadronization effects have to 
be included to further improve the theoretical description of the spectrum in the 
end-point regions. In this respect, the NLO calculation is vital to disentangle per- 
turbative effects from non-perturbative physics and relativistic corrections, and thus 
test the hadronization models that have been suggested in the literature. Finally, a 
more reliable NLO extraction of ag from radiative T decays is possible by restricting 
the analysis to the perturbatively calculable part of the energy spectrum. 
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